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ABSTRACT 

The Skyrme effective field theory is tested by evaluating nucleon ground state 
matrix elements of the correlation functions for two flavor density operators and 
two pseudoscalar density operators in the Skyrme model and comparing them with 
results in quenched lattice QCD. The possiblility of using quenched lattice QCD to 
study higher-order terms in effective field theory is also discussed. 
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I. INTRODUCTION 



The Skyrme model*^^"'^'' provides an appealing but incomplete step towards the goal of 
constructing a systematic, quantitative effective field theory for baryons. In the large N limit, 
QCD becomes equivalent to an effective theory of mesons, and baryons emerge naturally as 
solitons of this theory. The simplest Lagrangian for the effective theory, obtained by Skyrme 
by adding a non-minimal term to the non- linear sigma model to stabilize these solitons^, 
provides a remarkably successful schematic model of the nucleon. However, despite the large 
body of research studying specific corrections and extensions since Adkins, Nappi and Witten^ 
revived interest in the Skyrme model a decade ago by proposing it as a serious model of the 
nucleon, the effective field theory has never been developed and tested in a controlled way. 

Hence, the motivation for this present work is to use quenched lattice QCD as a frame- 
work for the controlled and systematic study of effective field theory. In contrast to studies 
of physical hadrons, for which the omission of fermion loops in quenched QCD represents an 
uncontrolled approximation, the fact that fermion loops are negligible in the large N limit 
implies that there is a meaningful effective meson theory for quenched QCD. Hence, lattice so- 
lutions of quenched QCD provide an ideal laboratory to systematically test every aspect of the 
approximations introduced in effective field theory. Relative to comparisons with experiment, 
the lattice QCD laboratory provides the opportunity both to explore the approximations 
as a function of fundamental parameters, such as the quark mass, and to study quantities 
which are not readily amenable to experiment, such as two-body correlation functions. Since 
quenched lattice QCD reproduces masses and hadron properties which are close to experiment 
for the physical quark masses, we expect the conclusions concerning effective field theory to 
be relevant to physical hadrons. 

To clarify the aspects of the effective field theory formulation we would ultimately like 
to explore, it is useful to consider the functional integral for full QCD. Meson fields may be 
introduced by writing delta-functions setting each meson field equal to appropriate bilinear 
quark fields and integrating over the meson fields. The relevant effective action is then defined 
by the functional integral over quark and gluon fields with the delta-function constraints, and 
this action is to be integrated over all the meson fields, corresponding to the calculation 
of all quantum loop corrections to the stationary-phase approximation. Similarly, effective 
operators are defined by the ratio of the functional integral over quark and gluon fields of 
the action multiplied by the original quark and gluon operator and delta-function constraints 
divided by the corresponding integral without the operator. In this framework, one could 
systematically study the effects of truncating to various numbers of meson fields, of truncating 
the terms in a derivative expansion of the effective action, of truncating at various levels in 
the loop expansion of quantum corrections, and of retaining various terms in the expressions 
for effective operators. In this language, the simplest version of the Skyrme model is obtained 
by retaining only the pion field, arguing by symmetry and simplicity that only the non-linear 
sigma model kinetic energy and Skyrme terms need be retained in the action, assuming that 
N is large enough that quantum corrections to the properly projected degenerate classical 
solutions may be neglected, and arguing that the relevant effective operators are uniquely 
specified by symmetry considerations. In principle, each of the four classes of approximation 
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may be tested quantitatively by lattice calculations, affording the opportunity to explore 
systematically all of the approximations underlying the effective theory. 

In this present work, we take the first step in such a program by comparing correlation 
functions in the nucleon calculated on the lattice with corresponding correlation functions 
calculated in the simplest version of the Skyrme model. For each value of the quark mass, 
we take the the lattice observables as defining a model system, which we approximate by an 
effective Skyrme Lagrangian. As in the work of Adkins and Nappi^, we use the values of the 
pion, nucleon, and delta masses to specify the three parameters in the Skyrme Lagrangian, 
which then predict all other observables. Since we wish to study the validity of the effective 
theory for hadron structure, we focus our attention on two correlation functions we have 
calculated previously which explore the spatial distribution of the nucleon: the density-density 
correlation function and the pseudoscalar correlation function. Naively, since the original 
Skryme model with parameters determined from hadron masses ^ makes a 40% error in Ftt, 
we expect errors up to this order of magnitude in other observables. In addition we expect 
the effective theory to improve as the pion becomes lighter. To the extent to which these 
expectations are borne out and the effective theory provides a reasonable first approximation 
to the lattice results, we believe we have established a useful framework for quantitative 
investigation of each of the corrections discussed above. 

The outline of this paper is as follows. In section II, we describe the operators and cor- 
relation functions which are evaluated on the lattice and our identification of corresponding 
operators in the Skyrme model. Using the standard spin-isospin projected hedgehog solu- 
tion, analytic expressions for the correlation functions are presented and the details of the 
derivations are given in the Appendix. The results for the density-density and pseudoscalar 
correlation functions in the Skyrme model are compared with corresponding lattice results in 
section III and the conclusions are discussed in the last section. 

II. Operators and Correlation Functions in Lattice QCD and in the Skyrme 
Model 

Since the QCD and the Skyrme Lagrangians are defined in terms of different degrees of 
freedom, in the absence of a systematic derivation of effective operators as discussed in the 
introduction, the determination of the operators in the Skyrme model to be compared with 
lattice operators might appear ambiguous. However, at the level of truncation considered here, 
we will show that the appropriate operators may be determined from symmetry considerations. 

The Skyrme Lagrangian has been constructed such that it shares those QCD symmetries 
believed to be relevant to low-energy phenomenology: the SU{2)l (S> SU{2)fi fiavor symme- 
try, with its SU{2)v subgroup that corresponds to isospin conservation, and a topologically 
conserved current that corresponds to baryon number conservation. Since each operator of 
interest to us in QCD is induced by a symmetry operation on the QCD Lagrangian, we de- 
fine the corresponding operator in the effective theory as the operator induced by the same 
symmetry operation on the Skyrme Lagrangian. 

density operators 

The two fiavor density operators we measure in lattice QCD are Pu{r) =:u{r)'y^u{r): 
and Pd(r) =:d{r)^^d{r):. These operators may be written as linear combinations of the 
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time-components of the conserved baryon current and conserved isospin current F^'" 
(o = 1, 2, 3), which correspond in the standard way to the isoscalar and isospin symmetries 
of QCD: 

h(^) = lB\r)±^-V''%r) . (2.1) 

Our normalization conventions are such that J d'^rV^'"'{r) = 2Ia (a**^ component of the 
isospin operator), J d^r B^{f) = B (baryon number operator) and J d'^r pu.di'f') = N^^d (flavor 
number operator). Since isospin is unbroken in QCD with equal quark masses and in the 
Skyrme model, both currents are exactly conserved. Hence they do not get renormalized and 
the normalization in Eq. (2.1) has an absolute meaning. Although it is possible to write 
conserved currents on the lattice, the local currents we have used are not conserved and 
have normalization factors which differ from unity by the order of 10%. As will be shown 
subsequently, this is not a problem in practice since the isoscalar current strongly dominates 
the isovector current. Therefore a 10% error in the relative normalization is negligible and 
the overall density is normalized such that it integrates to unity. 

We then may write the correlation between the up and down quark density as 

{pu{r)Pd{r')) = I {B\r)B'{r')) - \ (v'^\f)V'^'{f')) 

c^ss. 9 ^^0(^)^ l^^o^^.^\^ _ 1 ^^0,3^^)^ ^t>°'^(f')) (2.2) 

= -^B\f)B\f') - lF0'3^f)V0'3(f') . 

While the first equality is valid in general, the second one is only valid if we consider the 
classical solution and disregard the quantum fluctuations^. This is the approximation in 
which we solve the Skyrme model. 

We therefore compare the density-density correlation function evaluated in the nucleon 
ground state on the lattice 

/ \ ( \ - f d^rd^r'5{y- \r-f\) , 

{poPo)Liy) = J ^^^p {pu{r)Pd{r)) , (2.3) 

with the following function calculated in the Skyrme model 

{popq)s ^y) = i] ^ ) 

_1 j d^T d^r' 5{y - |r - f'\) 



(2.4) 



Here Bi^{r) and V'^'^{r) are the standard currents of the Skyrme Lagrangian: Bf^{f) is the 
topological baryon current^, and V^^'°' are the Noether currents corresponding to isospin 
symmetry. These currents are evaluated for the hedgehog solution and normalized such that 
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J d^rB^{r) = 1 and J d^rV^'^{r) = 2/3. The explicit forms of these currents, using the 
hedgehog solution to the classical equation of motion, U{r) = e*'""^^*^'"^ , are: 



A(f) 
A 



{eF^y sin'^ F{r) dF{f) 
27r^ df 

,A(r) 
A 



2In 



3/ • rflffl - la 

2 



sin2F(r) <^ 1 + 4 



d^rA{r) 



=3 1? 

^ J. TT 



df 



-A 



sin^ F(r) 



(2.5) 

(2.6a) 
(2.6b) 
(2.6c) 



The definition of A agrees with Ref. [2] and for convenenience we have defined the dimensionless 
constants A and f = eF^r. The vector fji is the unit vector in the f direction rotated by a 
time-dependent angle, where this rotation is introduced to project the hedgehog onto states 
of definite spin and isospin. Since we integrate functions that depend only on the relative 
direction between f and f' over the angles dO(f) and dQ{f'), we can always rotate back to 
f, f' and we never need to know the rotation explicitly. Details on how these currents and 
F{f) are derived are given in the Appendix. As mentioned previously, lattice measurements 
of the pion, nucleon and delta masses are used to determine the numerical values of the three 
parameters that enter in the equation for F{f): the pion mass ttt-tt, the pion decay constant 
and rho-pion coupling constant e. 
As shown in the Appendix, the final formula we obtain for nucleon states, once angular 
integration has been performed is: 



1 



(eF,) 



3 (PoPo)skyrme 



37r 1 



9 1 ^^sin^ F(r) dF{f) 
f df 



4A2 yji^y 



siT?F{f') dF{f') 

(Jjf 7~~ 



\y-r\ 



dffA{f) / df'f'A{f') 



y-r\ 



1 + 



df 



2rr' 



(2.7) 



where y = eF^^y and we used {Is)" 



I for nucleon states. 



Pseudoscalar Operators 

The other operators we measured on the lattice are the pseudoscalar density operators 
/5^(r) = •u(r)7^'u(r) and p^(f) = d{f)'y^d(f). For establishing the correspondence between 
operators in QCD and in the Skyrme model, it is useful to note that these operators are 
proportional to global chiral variations of the QCD Lagrangian containing a mass term that 
explicitly breaks the symmetry. For notational simplicity, for a chiral transformation i/j — > 
e*x75^ we will denote the chiral variation ^ by ^c- Then, disregarding for a moment the 
fiavor content, a chiral variation of the QCD Lagrangian yields ScJO,qct> ~ 2mq'4}^^'^ ■ Using 
the fact that the ratio mq/m^ is finite in the chiral limit (proportional to F^/ < ipip >), an 
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equivalent form for p^(r) written purely in terms of the Lagrangian and pion mass which has 
a well defined chiral limit is: 

Zniq 

As for the flavor densities, we may write the chiral variation as a linear combination of isoscalar 
and isovector terms. If we denote the variation of the QCD Lagrangian under an isoscalar 
chiral rotation by ScsJ^qcb and the variation under an isovector chiral rotation by Scv J^qcb 
we may write: 

pl,d = ciu,d'm~'^5csJ^QCD + bu,d'm~^5cv J^QCD ^ au,d'm~'^d^S'^ + hu,dm~'^d^V^'^ , (2.9) 

where ^'5 and V^""" {a = 1,2,3) are the corresponding Noether currents which are conserved 
in the chiral limit, and the a's and 6's are known constants which are finite in the chiral 
limit. We may then define the analogous operators in the Skyrme model to within the 
ambiguity of the choice of the mass term, for which we make the standard choice of C — 
|m^F^(TrC/ — 2). Since the two operators on the right-hand side of Eq. (2.9) are not 
conserved, in general one would have to calculate the renormalized operators to determine 
the correct linear combination. However, for the Skyrme Lagrangian ScsJ^Skyrme = by 
definition. That is, the variation is outside the SU{2) manifold where the theory is defined 
and there is no anomalous contribution as for the baryon number. Hence, in the Skyrme 
model both fiavor operators have only the isovector contribution 

Physically, this result corresponds to the fact that the the Skyrme model describes the 
— correlation as a pion-pion correlation, and from PCAC we recognize m~^dnV^'^ as a 
pion interpolating field. There is no contribution from pseudo-scalar isoscalar particles since 
the ?7 and ri' are not included in the model. One might argue whether the model makes a 
good approximation in keeping only the light pion degrees of freedom, but once we are given 
this model Lagrangian, the comparison with QCD is unambiguous. We clearly expect this 
approximation to improve as we lower the pion mass. 

Only the mass term £' = im^F^(Tr[/ — 2) in the Skyrme Lagrangian contributes to 
the chiral variation ScvJ^Skyvme' 

m-'^bcvC! = -^F^TviiUrs) = ^F^ sin F{r)rR ■ z (2.11a) 

where again fa is the rotated vector and, once we change variables from r to tr, we obtain 

Puir) ^ Pd{r) ^ F^ sin F{r) cose . (2.11b) 

Note that the operator has a well defined chiral limit in the Skyrme model. Hence, finally we 
compare the correlation function evaluated in the nucleon ground state on the lattice 

^^-1 '^'''''^%:^'-''^^ mm')) , (2.12) 
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with the following Skyrme model correlation function: 



(P5P5)s {y) 




d^rSr' 5 (y - \r - r'\) 



sin F{f) sin F(r') cos 9 cos 9 



(2.13) 



dr' sin F(f') [r^ + r'^ - y^] 



The overall normalization cannot be defined unambiguously as explained before, and we can 
only compare the shapes of the two functions. Details of the angular integration are given in 
the Appendix. 

To explore the ambiguity in the pseudoscalar current induced by the choice of the mass 
term, we have considered two alternatives to the standard result in Eq. (2.11) and compare 
the resulting correlation functions in Figure 1. Since we make the comparisons in the chiral 
limit, = 0, it is clear that each of the currents is consistent with the Lagrangian and simply 
corresponds to a different arbitrary choice of non-linear terms in the pion field. The standard 
result from Eq. (2.11) is shown by the solid curve and labeled by sin(F). Variation of the 
mass term C = ^m^F^ (Trt/^ — 2) yields a pseudoscalar density proportional to ^ sin(2F) 
and the resulting correlation function is shown by the dot-dashed line. In general, any term 
or combination of terms of the form Cj^ = gp-m^F^ Tr (U'' — 1) gives the expected asymptotic 
pion mass term and provides a possible alternative. The curve labeled by F is obtained by 
replacing sin F(f) in Eq. (2.11) by its asymptotic form at large r, F(f). Note that this result 
may be viewed as being generated by the Lagrangian term C — ^fe^i(— 1)'^>C^. This simple 
comparison shows explicitly that we can get very different results for correlation functions 
at short and medium range for different choices of effective operators which are equivalent 
at large distance. This clearly motivates future efforts to derive or calculate on the lattice 
higher-order terms in the expansion of effective operators. 

In a previous comparison^ of the lattice results and the bag model, it was natural to equate 
the quark field operators of the two theories. In the present context, one may ask whether the 
results would have been different if we had used operators induced by symmetries. It is readily 
seen that whereas we would have obtained the same fiavor density operators, the pseudoscalar 
density operators would have been somewhat different. The main point is that in contrast to 
the QCD and Skyrme Lagrangians, the bag Lagrangian in the basic form used in ref. [7] is not 
invariant under chiral rotation, even in the limit of zero quark mass. In the chiral limit, the 
chiral variation receives non-vanishing contributions from the bag surface. This contribution 
is equal in form to the QCD term, but it is multiplied by a delta function at the surface. 
Hence, in addition to the double volume integral of the relevant quark wave functions, we 
would also have had a double surface integral with a weight that would have become dominant 
in the chiral limit. The shape of this latter contribution is significantly different, even though 
it retains some of the correct global features such as having zero volume integral^. 
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III. RESULTS 



We use quenched lattice QCD calculations at three different values of the bare quark mass 
to define three test systems which may then be approximated by the Skyrme model. The 
lattice calculations are carried out on a 16^ x 24 lattice at /? = 5.7 using bag model sources and 
are described in detail in Refs. [7-10]. We have chosen to define the energy scale of the lattice 
calculations using the string tension, which gives an inverse lattice spacing = IGeV and 
a — 0.2fm. Although defining the scale using the extrapolated proton mass yields a value of 
a 15% smaller, the quantitative comparison with the Skyrme model is unaffected since when 
the Skyrme parameters are determined from the lattice masses, the correlation functions scale 
in a. 

The equation for the soliton, Eq. (A. 5), depends only on the dimensionless parameter 
P = and its solution is a function of the dimensionless variable f = cFt^t. Therefore, the 
correlation functions we calculate depend only on two parameters, m-j^ and the scale factor 
(cFtt)"^. However, since we link the choice of eFj^ to the proton and delta masses through 
Eq. (A. 3), we actually determine e and F^^ separately. The four sets of Skyrme parameters 
we use are shown in Table I, where the masses in bold face have been used as input and the 
other parameters are predicted by the model. Since the nucleon and delta energies depend 
nonlinearly on , the parameters are determined numerically. A rough idea of the sensitivity 

of the nucleon or delta mass to eF^ is given by the fact that at m-,^ = 0, ^ ~ 2^^^ 

density-density correlation function 

As a prelude to considering the Skyrme model density-density correlation function's de- 
pendence on the quark mass and how it compares with the lattice results, it is useful to note 
the extent to which it is dominated by the isoscalar baryon number density contribution. The 
isovector and isoscalar contributions, as well as the weighted average | < pp >5 — | < pp >v 
are shown in Figure 2 for the range of pion masses we consider. In all cases, over the spatial 
range shown in the figure relevant to this work, one observes that the isoscalar and isovec- 
tor contributions are sufficiently similar that | of their difference is very small compared to 
the isoscalar contribution and the full correlation function is very well approximated by the 
isoscalar correlation function. This result has two important implications. We may regard 
the physics as being dominated by the topological baryon number density and, as claimed 
earlier, we are justified in neglecting differences in the renormalization factors for the isovec- 
tor and isoscalar density operators on the lattice. The only place the isovector and isoscalar 
correlation functions differ significantly is in the extreme surface, where the isoscalar corre- 
lation decays asymptotically as e~^"^'^^, whereas the isovector correlation function decays as 
g-2m^r ^pj^jg behavior refiects the leading contributions of two or three pions respectively in 
each channel and gives rise to the fact that at very large distances, the correlation function 
changes sign. Similarly, we also expect the correlation to be negative at sufficiently large 
distance on the lattice. 

The change in the shape and spatial extent of the correlation function with quark mass 
is shown in Figure 3. For convenience, throughout this section we will refer to pion mass 
dependence interchangeably with quark mass dependence. It is implicit that when we consider 
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a change in the pion mass, we also consider the correlated changes in and e induced when 
a change in the bare quark mass produces this change in the pion mass. 

A striking feature of Fig. 3 is the fact that the mass dependence is highly non-linear. 
Indeed, the change in the correlation function caused by reducing m.j^ by 173 MeV from 515 
to 342 MeV is quite small compared to the change caused by reducing uIt^ by another 205 
MeV to 137 MeV. A related indication of how far these masses are from the chiral limit is 
given by the values shown in Table I for the parameter P = governing the hedgehog 
solution. 

The radial extent of the correlation function depends both on the scale factor (cFtt)"^, 
which decreases as m^r decreases, and on the exponential tail, which grows as 171^^ decreases. 
These two effects go in opposite directions for the cases we consider. It turns out that the 
scale factor dominates and the net result is shrinking of both the isoscalar and isovector 
contributions with decreasing m-j^. This shrinkage with decreasing is shown both by 
the increase of the normalized correlation function at zero separation and by the decrease 
in isoscalar and isovector r.m.s. radii in Table I. (Note that for values of 111^^ lighter than 
the physical mass, the isovector radius begins to increase and ultimately diverges in the 
chiral limit because of chiral logs.) The trend is opposite in the lattice results, where the 
normalized correlation function at zero separation clearly decreases with m^r. The second 
moments of the lattice results have sufficiently large errors because of the weighting of the 
poorly determined tails that their trend with m-j^ is not statistically determined. Also note 
that these second moments of correlation functions are not rms radii of density distributions 
because of the presence of polarization contributions^^. 

Lattice measurements of the density-density correlation function for each of the three val- 
ues of the bare quark mass are compared with the corresponding Skyrme model correlation 
functions in Figure 4. Because statistical errors in the correlation functions at large radius 
induce significant errors in the overall normalization, we have normalized all correlation func- 
tions to unity at the origin to facilitate comparison. One observes that as the quark mass 
and thus the pion mass decreases, the agreement improves substantially. Since we are using 
sufficiently large quark masses that we should not expect the Skyrme model to be an accurate 
approximation, it is encouraging that the agreement improves so conspicuously as the quark 
mass decreases. It will therefore be of interest to carry out more detailed comparisons for 
lighter quark masses, to verify that this agreement is maintained as the correlation function 
undergoes the significant changes shown in Fig. 3. 

pseudoscalar correlation function 

In contrast to the density-density correlation function, where both isoscalar and isovector 
terms contribute, the Skyrme model has only an isovector contribution to the pseudoscalar 
correlation function. This correlation function decays exponentially in the pion mass at large 
r and, on physical grounds, it should be the dominant contribution in the chiral limit at large 
distance. In contrast, the lattice calculation has contributions from both the isoscalar and 
isovector terms. At the large quark masses we consider in this work, the isovector tt and the 
isoscalar rj have comparable masses and may be expected to yield comparable contributions 
with an overall coefficient that favors the isoscalar one. Hence, we expect the best agreement 
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between the Skyrme and lattice results to occur near the chiral limit, where the isovector 
contribution becomes dominant. 

The dependence of the spatial distribution of the pseudoscalar correlation function on the 
quark mass is shown in Figure 5. Although in principle one has the same competition between 
the shrinking scale parameter and growing pion tail as in the density-density correlation 
function, we observe that the pseudoscalar correlation function has extremely mild mass 
dependence in the range 0-2 fm. The only qualitative feature which is not shown in this 
plot is the fact that in the Skyrme model, the volume integral of this correlation function 
is zero, as it is in any independent particle approximation^. Thus, the function has a node 
and an exterior negative region of equal volume. When the correlation function is plotted at 
larger distances, one observes that as the pion mass becomes lower, the tail extends to larger 
distances and the node moves outward. 

The comparison of the Skyrme model pseudoscalar correlation functions with the lattice 
results is shown in Figure 6. On the lattice we observe that the overall size is very different, 
as expected because of the large isoscalar contribution at the quark masses we have used. 
However, it is significant that as the pion mass decreases, the lattice result becomes much 
closer to the Skyrme solution. The outward movement of the node is particularly evident in 
Figure 6. The stability of the Skyrme correlation function when varying the mass compared 
to the large change of the lattice result suggests that a significant part of change in the lattice 
result arises from the diminishing contribution of the isoscalar term. Again it would be very 
interesting to follow the behavior to lower quark masses. As a technical point, we note that 
because of the node at large distance and the slow spatial decay at small pion mass, the 
pseudoscalar correlation function is more sensitive to finite size effects on the lattice, and is 
more difficult to extrapolate to the physical pion mass than the density-density correlation 
function. 

IV. SUMMARY AND CONCLUSIONS 

In this work, we have derived the analytic results necessary to compare the density-density 
and pseudoscalar correlation functions calculated in lattice QCD with the Skyrme model and 
taken the first step in using quenched lattice QCD as a laboratory to systematically explore 
the Skyrme effective field theory. The principal numerical results are presented in Figures 4 
and 6. We observe in Figure 4 that already for a quark mass of 40 MeV, corresponding to a 
pion mass of 340 MeV, the Skyrme approximation to the density-density correlation function 
appears to be in good agreement with the lattice result. In contrast, because of substantial 
isoscalar contributions in quenched QCD for a pion mass of 340 MeV, which are excluded 
by construction from the Skyrme model, the Skyrme pseudoscalar correlation function is not 
quantitatively correct, although the trend with decreasing pion mass is encouraging. 

These exploratory results suggest several promising directions for future research. At 
the most pedestrian level, it will clearly be worthwhile to extend the comparison to lighter 
quark masses. More interesting and potentially extremely useful, is the opportunity to use 
lattice techniques to calculate numerically the effective action and effective operators and 
thus evaluate corrections to the Skyrme approximation. As described in the introduction, 
this provides a consistent framework to systematically explore the quantitative effects of 
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truncation in meson fields, of truncation of the effective action, of truncation of the expansion 
of effective operators, and of truncation of the quantum loop corrections. Thus, rather than 
using the lattice to calculate hadron observables directly, one can instead use it to calculate 
from first principles the relevant parameters of an effective theory which could then be applied 
much more generally than one could hope to apply the full apparatus of lattice QCD. 
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APPENDIX 



In this Appendix we review the main features of the Skyrme model, we give an exphcit 
derivation of the time component of the vector current without averaging over angular vari- 
ables (which is also useful for other applications) and we perform the angular integration 
when the two currents are kept at fixed relative distance. As in the rest of the paper we have 
tried to keep notation consistent with Refs. [2,3], from which Eqs. (A.l) through (A. 7) are 
taken. 

The Lagrangian of the model is (Ref . [3] ) 
jC=j^F^tT[d,Ud^U^+^tT[[{d,U)U^,{d,U)U^]']+lmlF^itrU-2) . (A.l) 
If we substitute the spinning hedgehog ansatz 



U = A{t)Uo{f)A\t) , A{t)eSU{2) 

- (A. 2a) 

C/o(r) = e^-^(^)^"" 

into (A.l), we obtain after quantization the following expression for the energy (Ref. [2]): 

(A.3) 

where J is the spin (equal to isospin in this model) of the particle 



e 



M+ ^2J(2J + 2) 
8A 



M = I d'f i - 



— 2 

dr 



sin^F 



1 

+ 2 



dF\^ sin^F 



2(-pl + 

ar 



sin^ F 1 ^2 / 



(A.4) 

A was defined in Eq. (2.6c), r = eFj^r and /3 = F{r) is chosen so that it minimizes M 

with the boundary conditions F(0) = tt and F(oo) = to ensure that the baryon number 
(integral of Eq. (2.5)) is equal to 1; the resulting equation for F(r) is 



-r^ + 2sm^F — r- + 2smFcosF — + -r — 
4 J df^ \dr J 2 dr 



- 2sinFcosF - + 



1 sin^ F 



(A.5) 



- ^/^VsinF = 
4 



The general form of the time component of the Noether current associated with the V — A 
transformation 5U = iTaU of the Lagrangian (A.l) is (Ref. [2]): 



4- A = ^iF> {{d'U) t/tr4 + ^tr {[(a.t/)t/t,r,] [{d'U) U\ {d^U)U^]] 
= liF^ tr { (d'U) U^Ta] - g^E, tr { [{diU) U^Ta] [{d'U) U\ {d^U) U^] } 



(A.6) 
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By using Tr {[A, B] [C, D]} = 2 Tr {AD} Tr {BC} -2Tr {AC} Tr {BD}, valid for A, B, C, D 
belonging to the SU{2) Lie algebra and the fact that d^U^ = d^U~^ = —U~^ i^nU) = 
—U^ (dfj^U) C/^, we may rewrite (A. 6) as 



7-0,0 _ ■^-^■^ 



tr {{d'U) U^Ta} - ^tr {{d'U) UWa}^^tr[[{^,U)U^'} 
+ ^S,tr {idiU)U^Ta} tr {{diU)U^ idoU)U^ 

?2 



^tr {{doU) U^Ta} + ^tT U^Ta}^itT {{diU) {diU^)} 

-^Sitr {{diU)U^Ta} tr {{diU) {doU^)} . 



(A.7) 



4e2 

Before calculating the four traces we need, let us derive a few useful formulae. Equation 
(A. 2a) can be written as 



Uo (r) = cos F + if • f sin F 
C/(f , t) = A{t)UoA^ (t) = cos F + if ■ tr sin F 



(A.2b) 



where tr is the vector f rotated according to the law tr - r = Ar • tA'^. Using the fact that 
A is unitary, we obtain 



^ = u^ AuqA^ + auqA^ = Aa^auqA^ - auqA^Aa^ = \Aa^, u 

at I 
The spatial derivatives are 

dUo ^ ^ dF sinF . . 

— — = diUo = ri — (- smF + zcos Fr • r) + t [Ti - rir ■ r) 

oxi dr r 

dU ^o^^ , dF . . ^ . ^, .sinF , ^ ^ ^, 

— ^ = [/ = "^Ri^ {-smF + icosFrR ■ r) + ^ [n - T^RiXR ■ r) 



(A.8) 



(A.9) 



We also need to show that 



or equivalently that 



s,[a,(...)] m...)] = ^i [df (. . .)af (. . .)] 



(A.lOa) 



(A.lOb) 



This is because it is more convenient to take derivatives of U with respect to tr (see Eq. (A.9)). 
Since 

we may write 



Jjk = Js^tr {n {A^TjA)} tr {n (^Vfe^)} 



(A.lOc) 



where we used the fact that 



1 tr {Ml M2} = ^ tr {Ml} ^ tr {M2} + S i tr {r^Mi} ^ tr {r^Ma} , Mi, M2 G 5t/(2) . 

(A.ll) 

This is the component form of the scalar product in SU{2). Now we are ready to calculate 
the traces in Eq. (A. 7). The first one gives 

tr {{doU) U^a} = tr I [iAt, C/] U^a} 

= tr I AA\ i sin Ff_R • fj (cos F — z sin Ffn • f ) Taj 
= zsinFcosFtr | |^^^''',ffl ■ rj r^j + sin^Ftr | |^lA't',fii • fj fi? • TTaj 
= isin(2F) tr ji^Vfl • fra} + 2sin2 Ftr ji^Va} 
- sin^ F tr |i^V^-f} tr {fii-TTa} , 

(A.12) 

where we used ir {[A,B]CD} = tr {^D} tr {SC} - tr {^C} tr {SD} and tr{^SC} = 
— tr {BAG} for A, B, C, D, belonging to the SU (2) Lie algebra to go from the second last to 
the last line. The second trace is 



:tr {(diU) {diU^)} = S^tr {(a^t/o) (SiC/^) } 

^ ( \ ^ dF ^^ sinF . . 

tr i Ti—— (— sm F + I cos Fr • r) + t [Ti — r^r • r) 

dr r 



X 



, dF ^\ smF , „ _x 

rj— — (— smF — t cos Fr ■ r) — ^ (r^ — r^r ■ r) 

ctr r 



/dFV sin^^ 



= tr < — + 



f - 2 + 1 n = 2 



fdFy ^sin^F 
\dr J r^ 



\dr J 

(A.13) 

The third trace is, using the fact that in the last term of Eq. (A. 7) we shall take advantage 
of property (A. 10a) to change di into 9/^, 



tr {(afC/) U^Ta}= tr 

X [cos F — z sin Ff_R • f ] | 



\ dF . . ^ . .sinF, ^ ^ ^' 

rm — (- sm F + z cos Ftr ■ r) + i (rj - rmrR ■ r) 

dr r 



(A.14) 



sin -Z^ 

= TRi tr {...} + z Tr {[cos F - z sin Ftr ■ r] Tan} ■ 



We do not write out the term proportional to tri, because it will not contribute when multi- 
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plied times the fourth trace, which is 



tr {{doU) (af c/t)} = tr I \Aa-\ U 

= tr jiA-^ 
sin^F 



afC/-i} = tr [U,d^U-^]] 



^ .sinF 
isinFru ■ r, -i 



tr [aA-^ [tr ■ f, n]] = tr [AA'^tr ■ rn} . 

(A.15) 

Now we multiply Eq. (A. 14) times Eq. (A.15) and sum over i to get the last term of Eq. (A. 7). 
Note that the first term in Eq. (A. 14) when multiplied times Eq. (A.15) gives a result propor- 
tional to tr = 0. To combine the trace in the second term 
with the one in Eq. (A.15) we may use Eq. (A. 11): 



S,tr {{d,U) U'Wa} tr {{d,U) {doU^)} 

= Sitr {(aft/) u^Ta} tr {(aft/) (aot/t)} 

sin"^Fl f ; . + . w 1 1 rr . -, 

--tr <^AA^rR - TTij - tr {[cos F - ism FrR ■ t\ Tan} 

tr l^iylVit • [cos F — i sin Ftr • f ] | 



z8- 



= i8- 



sin-^F (\ 



(A.16) 



tr l^dA^ffl • ^} ^ tr {-isva-Fr-R ■ TTa}^ 

2^^^2— [isin2Ftr ji^Vfl • fra} + 2 sin^ Ftr jiAVo} 
- sin^ F tr \aA^tr • f } tr (fi? • rrj] . 



The expression in brackets is equal to tr {(aot/) t/Va} (Eq. (A. 12)), which appears in both 
the other two terms in Eq. (A. 7). We then insert Eq. (A. 12), Eq. (A. 13) and Eq. (A. 14) into 
Eq. (A. 7) and obtain 



0,a 
V-A 



iFl 



sm 



Ftr {iAV„} + ^ sin2Ftr \aA^r -tt^- ffiaSin^ Ftr \aA^r ■ f } 



X 



1 + 



(eF^)2 \ \dr ) 



dFV sin^F 
' + T- 



(A.17) 

As in Ref. [2], one can convince oneself that Jy+A obtained from Eq. (A. 7) by exchanging 
U <-> U~^, i.e. "Tr <-> —"Tr. Then Jy will contain terms even in tr and Ja terms odd in tr, 
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as expected from parity considerations: 

iF? 



jOjU TO,a . jOjU 



sinV 



1 + 



X [tr -fflatr |i^^ffl-rj 

= 3iF^A{f) [tr { iA" Va } -fiia tr ji^Vji ■ r} 



sin^F 



(ir 



+ 



(A.18) 



p2A(r) 



{2/a 



3/ • rufiia - I. 



where we used the definitions of A and A, Eqs. (2.6b) and (2.6c). We also used the fact that 
AA-^ = -2^/-r with / the isospin operator, as may be derived from A = ao + ia-T, tt = 4Xa 

and / = — I [aoTT — airo — a x tt] which can be found in Ref. [2]. For completeness, the time 
component of the axial current is 



p2 ( 4 

jO,a _ jO,a _ TO,a — -^^r ■ OF j 1 I 

-J A - Jv+A Jy-A - 4;^ sm^i^ < i + ^^^^^2 



(dFV sin^F 
\dr J 



{ixru)^ , (A.19) 



in agreement with Ref. [3]. 

We calculate only the correlations between currents in proton (or neutron), {N\VV\N), 
and not off-diagonal matrix elements, {N\VV\N'). Hence, we only need: 



diag 4''{r)4'{r') 



^ diag 



{I3-I ■ rrs) (is-I- r'r'3 



= F^4 9A(r)A(rO ^^^^^^ _ ^^^^ ^.^^ ^^^^^^ _ ^.^^ ^^^^^^ ^ ^.^^ ^^^^^^ 

^ ^4 9A(r)A(rO ( ^3)2 | ^ _ ^^^2 ^ _ ^^^2 ^/ ^ ^ ^ ^/ 1 ^ 

(A.20) 

where 7 is the angle between r and r' . 

The correlation we want is then (note that we changed variables of integration from f, r' 
to TR^f'ji and that (/a)^ = \ both in the proton and neutron): 



V{y) = J 



4 



dfdf 6{y — \r — r 



diag 4\r)4\r') 



9 1 



-1 rOO rOO p p 

^—^ j dff^j df'r''^K{r)K{r') j dcj) j dcj)' 



4 47ry2 (e^F^r 

X j d{cos9) J d{cos 9')S{y - \r -r\) [1 - cos^ 9 - cos^ 9' + cos ^ cos 9 cos 9'] 



Now we use 



6 (y — \r — r'\^ — 2y6 {y^ — (r — r^)^) = ^ ^ 6 | C0S7 — ^ ^ 



rr' 



2ff' 



(A.21) 
(A.22) 
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where y = eF^; we change variables from 0, ^, 0', 9' to 0, ^, 7, where ^» and 7 are the angles 
of f' with respect to r (the transformation is just a rotation and its Jacobian is one); and we 
use the following identities 

cos6'' = cos6'cos7 ± sin^sin7cos(0 ± $) , (A. 23a) 



/ 



d(j) cos 9' = 27r cos 9 cos 7 , 



/ 



(i(^cos^ 9' = 27r 



cos^ 9 cos^ 7+0 ^^^^ ^ ^^^^ 



(A.23b) 
(A.23c) 



Note that we need only (A. 23b) and (A. 23c) that can be derived from (A. 23a) without 
knowing the signs which depend on conventions. Then, 

9 1 [eF^Y 



4 [e^F^Xy 4Tvy 



' dfrA{f) / dr' f' A{r') {27: f d{cos9) / cZ(cos7) 
Jo J J 



X 5 ( cos 7 



X 



,^2 j y^^2 

2ff' 



1 - cos^ 9 - cos^ 9 cos^ 7 - - (1 - cos^ 6*) (1 - cos^ 7) + cos^ 9 cos^ 7 



1 



{eF^y 



X 



4 (e3F,A)2 y 

^2 _|_ ~/2 _ ~2\ 



poo poo 

/ dffk{r) / df'f'k{f') 
Jo Jo 



(l + z^) I dx{l-x^) 



2rf' J 2 

37r {eF^f 1 r ^~~M~\ H ~' K(~'\ 



X 



= 37r 



1 + 



2^' 



9 (2ff 



™/ 1^2 , ~/2 ~2 



r + r 



{eF^y 



y {e^F, 



— / dffk{f) / df'f'k{f') 



^2 _j_ ^/2 



(A.24) 

where in the last line we use the fact that the integrand is symmetric in r and r' to impose 
r > r' and multiply the result times two; with this condition 2rf' > \r^+r'^—y^\ r' > \r—y\ 
which in turn implies r > ^y. The other two angular integrals we need can be derived in the 
same way: 



/ 



y 



and 



/ 



47ry"^ 

dr dr 5{y — \r — r'\) 



drrF{r) / dr'r'F{r') 



(A.25) 



y-r\ 



Airy 

1 f dr dr 5{y — \r — r'\) 
3 

4-71 



F{r)F{r') cos 9 cos 9' 



/ 



3y 



47ry^ 

drrF{r) 



F{r)F{r') cos 7 



(A.26) 



2y 



dr' r' F{r') 



y-r\ 



J.2 _|_ _ y2 

2rr' 



Note that the hmit for vanishing y of each of the last three equations is finite^. 
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Figure Captions 

Fig. 1 Pseudoscalar correlation function in the Skyrme model for three alternative definitions 
of the pseudoscalar current. The solid curve denotes the standard current, Eq. (2.11b), 
induced by the usual choice of the mass term. The dashed and dot-dashed curves show 
two other alternatives discussed in the text. 

Fig. 2 Relative contribution of isoscalar and isovector contributions to the density-density cor- 
relation function in the Skyrme model. Dashed and dotted lines show the isoscalar and 
isovector contributions respectively at three representative values of the quark mass, and 
the weighted combination, Eq. (2.2), is shown by the solid curve. 

Fig. 3 Mass dependence of the density-density correlation function in the Skyrme model. The 
spatial distribution of the correlation function is shown for each of the four sets of pa- 
rameters in Table I. 

Fig. 4 Comparison of density-density correlation functions on the lattice and in the Skyrme 
model. The solid curves denote the Skyrme results for the parameters in Table I corre- 
sponding to each lattice quark mass. Lattice measurements of the correlation functions 
are shown by error bars joined by dashed lines to guide the eye, and all correlation 
functions are normalized to one at the origin. 

Fig. 5 Mass dependence of the pseudoscalar correlation function in the Skyrme model. The four 
curves show the spatial distribution of the correlation function for each of the four sets 
of parameters in Table I. 

Fig. 6 Comparison of pseudoscalar correlation functions on the lattice and in the Skyrme model. 
The solid curves denote the Skyrme results for the parameters in Table I corresponding 
to each lattice quark mass. Lattice measurements of the correlation functions are shown 
by error bars joined by dashed lines to guide the eye, and all correlation functions are 
normalized to one at the origin. 
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TABLE I 


Skyrme 


Lattice 


K 










0.167 


0.1639 


0.16 


ruq (MeV) 




40 


95 


175 


40 


95 


175 


(MeV) 


137 


342 


515 


694 


340(7) 


511(5) 


691(3) 


e 


4.83 


3.79 


3.23 


2.73 








(MeV) 


108.1 


60.32 


46.87 


34.87 








(eF^)-i (fm) 


0.383 


0.862 


1.321 


2.141 








/? = m^/{eF^) 


0.262 


1.497 


3.402 


7.290 








I IVic V 1 










Q1 ^(f\\ 


1 (^Q7(^ 1 ^ 




mA (MeV) 


1232 


1091 


1223 


1403 


1091(23)* 


1223(11) 


1403(10) 


rs (fm) 


0.690 


0.917 


0.970 


1.075 








rv (fm) 


1.048 


1.088 


1.132 


1.256 








f rs - \rv (fm) 


0.632 


0.894 


0.947 


1.050 


1.0(5) 


0.79(20) 


0.75(15) 



Table I. Parameters of the Skyrme model for lattice calculations corresponding 
to three different bare quark masses and for physical hadrons. The input hadron 
masses are shown in bold face font and all other quantities are derived from them. 
The isoscalar and isovector r.m.s. radii are denoted rs and ry respectively. The 
first column uses physical hadron masses and the next three columns use the lattice 
hadron masses given in the lattice data shown at the right. The lattice data are 
taken from Refs. [7-10], with the exception of uia at k = 0.167 denoted by an * 
which was taken from Ref. 12. For reference, we have also tabulated the values of 
the bare quark mass rUq = ^ — corresponding to the values of the hopping 
parameter k for each lattice. 
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12 3 



r (fm) 



Fiq„ 1 



1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 r 




= 694 MeV 




r (fm) 

Fiq„ 2 







1 1 1 1 1 1 1 1 1 1 1 1 


1 1 


1 1 1 1 1 


in^ - 


694 


MeV 


= 


515 


MeV : 


= 


342 


MeV 


- 

^^^^ t 


lor 


MeV 

- 






- 
- 


- — ^ ^ 


















1 1 1 1 1 1 1 1 1 1 1 1 


1 1 


1 • m-^^uTTTt*^ 







0.5 



1.5 



r (fm) 



Fiq= 3 




Fiq„ 4 



1.5 I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — r 




I I I I I I I I I I I I I I I I I I I I I 

0.5 1 1.5 2 

r (fm) 



Fiq„ 5 




Fiq= 6 



